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Radiatively inefficient accretion flows (RIAFs) in low-luminosity active galactic nuclei (LLAGNs) 
have been suggested as cosmic-ray and neutrino sources that may largely contribute to the observed 
diffuse neutrino intensity. We show that this scenario naturally predicts hadronic multi-TeV gamma- 
ray excesses around Galactic centers. The protons accelerated in the RIAF in Sagittarius A* (Sgr 
A*) escape and interact with dense molecular gas surrounding Sgr A*, which is known as the central 
molecular zone (CMZ), and produce gamma rays as well as neutrinos. Based on a theoretical model 
that is compatible with the IceGube data, we calculate gamma-ray spectra of the CMZ and find that 
the gamma rays with > 1 TeV may have already been detected with the High Energy Stereoscopic 
System, if Sgr A* was more active in the past than it is today as indicated by various observations. 

Our model predicts that neutrinos should come from the CMZ with a spectrum similar to the 
gamma-ray spectrum. We also show that such a gamma-ray excess is expected for some nearby 
galaxies hosting LLAGNs. 


I. INTRODUCTION 

Recently, the IceCube Collaboration reported the dis¬ 
covery of extraterrestrial neutrinos [iH^. The origin of 
the neutrinos is a matter of a debate (|4l4l9l| for reviews, 
see ©El). The data so far are compatible with an 
isotropic distribution, which suggests that neutrinos are 
of extragalactic origin. Diffuse gamma-ray data also sup¬ 
port this idea [1, [l^. However, the sources have not 
been identified because of poor angular information and 
statistics for the neutrinos. One way of improving this 
situation may be detection of counterparts through elec¬ 
tromagnetic waves. 

Low-luminosity active galactic nuclei (LLAGNs) are a 
candidate for the source of the neutrinos. The LLAGNs 
are expected to have radiatively inefficient accretion flows 
(RIAFs), which are realized when the mass accretion rate 
into the supermassive black hole (SMBH) is relatively 
small (M/MEdd ^ 0.01-0.1), where MEdd = TEdd/c^ is 
the Eddington accretion rate, and LEdd is the Eddington 
luminosity In the tenuous and turbulent plasma 

in the RIAEs, cosmic ray (CR) protons may be accel¬ 
erated via stochastic acceleration or magnetic reconnec¬ 
tion [l^ . These CR protons interact with other nucleons 
{pp interaction) and photons (pj interaction) in the flow 
and generate neutrinos. Although the production rate 
of the neutrinos per an LLAGN is not large compared 


with other more energetic sources such as quasistellar 
objects (QSOs), the abundance of LLAGNs can repro¬ 
duce observed neutrino flux on the Earth M- Even if it 
is difficult to resolve LLAGNs as point neutrino sources, 
gamma rays from them could be used to test this model. 
In particular, the gamma rays that are a byproduct of the 
pp interactions have energies comparable to those of the 
neutrinos, which means that the gamma-ray spectrum 
should reflect the neutrino spectrum unless the gamma 
rays are not absorbed. 


Sagittarius A* (Sgr A*) is the SMBH at the center of 
the Galaxy and it is known as an LLAGN. The current 
mass accretion rate of Sgr A* is very small and the accre¬ 
tion flow is thought to be a RIAF [23j. The current pro¬ 
duction rate of CR protons in the RIAF is expected to be 
small because of the small accretion rate [2^. Thus, the 
gamma-ray luminosity of the RIAF in the TeV band is 
also expected to be small because of inefficient pion pro¬ 
duction [l^. However, it has been indicated that Sgr A* 
was much more active in the past 25-2^. During those 
activities, a large amount of protons could have been ac¬ 
celerated and escaped from the RIAF. Moreover, observa¬ 
tions have revealed that there is a huge amount of molec¬ 
ular gas surrounding Sgr A*. This gas concentration is 
known as the central molecular zone (CMZ) with the size 
of i?cMZ "^100 pc and the mass of Mqmz 10^ Mq [29| . 
Strong turbulence and magnetic fields in the CMZ may 
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delay the diffusion of the CR protons that have plunged 
into the CMZ, and those protons may stay in the CMZ 
for a long time. In this paper, we calculate the diffusion 
of the protons in the CMZ that have accelerated and 
escaped from the RIAF in Sgr A*. We estimate gamma- 
ray and neutrino emissions created through pp interac¬ 
tions between the CR protons and protons in the CMZ. 
We show that TeV gamma rays from the CMZ around 
Sgr A* and nearby LLAGNs observed with the High En¬ 
ergy Stereoscopic System (HESS) may be produced by 
this mechanism. Previous studies have shown that the 
gamma-ray excess observed in those objects at ^1 TeV 
cannot be explained by one-zone leptonic models [31)1 - 132] . 
Our model can resolve this issue, although there are also 
studies trying to explain the HESS observations b y pp 
interactions in the CMZ by different approaches | 32 l 436 l |. 


II. COSMIC-RAY PROTON ACCELERATION IN 
RADIATIVELY INEFFICIENT ACCRETION 
FLOWS OF SGR A* 


In our model, protons are accelerated in the RIAF of 
Sgr A*. Since the acceleration is confined in a small re¬ 
gion on a scale of a few tens of the Schwarzschild radius 
of the SMBH, we consider the acceleration based on a 
one-zone model as previous studies [3113, Hi. Accord¬ 
ing to the model of Ref. 0, the typical energy of the 
accelerated protons is determined by the balance between 
the acceleration time of the protons in a RIAF (Gcc.r) 
and their escape time from the RIAF. The escape time 
is comparable to the diffusion time of the protons in the 
RIAF (tdiff.R)- Thus, the Lorentz factor corresponding 
to the typical energy is obtained by solving the equation 
of tacc,R = idiff,R and the result is 
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where nip is the proton mass, m is the normalized accre¬ 
tion rate m = M /M^dd, a is the alpha parameter of the 
accretion flow [39|; C is the ratio of the strength of tur¬ 
bulent fields to that of the nonturbulent fields, /3 is the 
plasma beta parameter. Race is the typical radius where 
particles are accelerated, and Rs is the Schwarzschild 
radius of the black hole 0 . For the parameters of 
the RIAF, we take a = 0.1, C = 0.05, /3 = 3, and 
Race = 10 Rs as fiducial parameters, because the neu¬ 
trino flux at cv. 10-100 TeV obtained with IceCube is 
well reproduced by them 0 . 

We assume that the luminosity of the protons accel¬ 
erated in the RIAF is Rp,tot = iJcrMc^, where i/cr is the 
parameter and we take r/cr = 0.015 as the fiducial value 
following Ref 0|. When only stochastic acceleration is 


effective, the production rate of protons in the momen¬ 
tum range p to p -I- dp is 

N{x)dx oc , (2) 


where x = p/peut, Ku is the Bessel function, and b = 
3/(2 — q) | 40 l |. The power-law index of turbulence re¬ 
sponsible for the acceleration is assumed to be g = 5/3 
(Kolmogorov type). The cutoff momentum is defined as 
= (2 - g)^/(^“'jVeq = Peq/ 27 , where Peq = Rp,eq/c 
[lOl koj l. We determine the normalization of Eq. m so 
that the total power of the protons is Rp,tot- 


III. DIFFUSION OF PROTONS IN THE CMZ 


Protons accelerated in Sgr A* leave the acceleration 
site (RIAF) and disperse into the interstellar space. Some 
of them would enter the CMZ surrounding Sgr A*. We 
solve a diffusion-convection equation for the CR protons 
in the CMZ. For the sake of simplicity, we solve a spher¬ 
ically symmetric equation: 
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where / = f{t,r,p) is the distribution function, r is the 
distance from the Galactic center, p is the momentum of 
particles, k is the diffusion coefficient, u is the velocity 
of the background gas, and Q is the source term for the 
particles (Sgr A*). We assume that m = 0, because we are 
interested in the CRs inside the CMZ, which is too heavy 
to be moved by possible outflows from Sgr A*. We do 
not consider CRs carried by the outflows without entering 
into the CMZ. We assume that the CMZ is uniform and 
its dense gas occupies at r < Rcmz- 

The actual CMZ has a disclike structure and does not 
entirely cover Sgr A* 0. Thus, we expect that most 
of the CR protons do not plunge into the CMZ, and we 
assume that only a fraction A of the protons accelerated 
in the RIAF are injected into the CMZ. Thus, the source 
term in Eq. (I3|) is written as / Anep^Qdp = XLp^tot = 
XpcrMc^. Since the size of the CMZ (Rcmz 100 pc) 
is much larger than that of the RIAF, we treat Q as a 
point source. 

We assume that the diffusion coefficient of CRs outside 
the RIAF is given by 
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where Ep is the particle energy. This coefficient is for 
the ordinary interstellar space in the Galactic disc M|. 
We only consider resonant scattering, which is valid at 
sufficiently low energies. Although the diffusion coeffi¬ 
cient in the CMZ is not known, we apply Eq. o to 
stronger magnetic field cases. If there are strong mag¬ 
netic fields (R ~ mC) in the CMZ [0, the coefficient is 
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much smaller than that in the intercloud space around 
the Galactic center {B ^ 10/rG) [i^. In fact, it has been 
indicated that the diffusion coefficient in molecular gas 
around supernova remnants is much smaller than the or¬ 
dinary value . From now on, we fix magnetic fields at 
B = 1 mG for r < i?cMZ and 10/rG for r > i?cMZ- We do 
not consider stochastic acceleration in the GMZ. This is 
because the diffusion coefficient we assumed in Eq. o is 
too large for effective particle acceleration. In fact, previ¬ 
ous studies have shown that the diffusion coefficient must 
be as small as that for the Bohm diffusion for particles 
to be accelerated up to ^ TeV 34|, [4^ . It is not certain 
whether such a small coefficient is realized by turbulence 
in and around the GMZ. We do not include cooling of 
the protons in Eq. ©, because the cooling time is larger 
than the diffusion time estimated based on Eq. o (see 
later). 

GR protons interact with protons in the GMZ. For a 
given distribution function /, we calculate the produc¬ 
tion rate of gamma-ray photons usin g th e code provided 
by 45j and the formula provided by [4^ for Ep <1 and 
Ep > \ TeV, respectively. The results are not sensitive 
to the boundary energy (1 TeV). We also calculate the 
neutrino production rate at the same time. We consider 
the attenuation of very high energy gamma rays by pair 
production on the Galactic interstellar radiation field us¬ 
ing the results shown in Fig. 3 of Ref. T7|. However, 
the following results are not much affected by the atten¬ 
uation. The energy density of interstellar radiation field 
(~ 10 eVcm“^) is much smaller than that of the assumed 
magnetic field (~ mG) [i^. Thus, the gamma-ray emis¬ 
sion via inverse Gompton scattering by secondary elec¬ 
trons can be ignored. 


IV. GAMMA RAYS FROM THE GMZ AROUND 
SGR A* 


Observations showed that the radius of the GMZ at 
the Galactic center is i?CMZ,obs = 200 pc, the thick¬ 
ness is Hcmz obs = 75 pc, and the mass is Mqmz = 
2x10^Mq Thus, the average density is pcMZ = 

McMz/(7ri?cMZ,obs-f^CMZ.obs) = 1.4X 10“2^gcm“3. Since 
Eq. ([3|) assumes spherical symmetry, we define an effec¬ 
tive radius 

flcMzs(|^^^)‘'’ = 130pc, (5) 

\4 7rpcMZ/ 


and we use this as the radius of the GMZ in 
the following calculations. The covering factor 
of the GMZ for r < i?CMZ,obs is /cmz = 
7r7?cMZ.obs77cMZ.obs/(47ri?^MZ.obs/3) = 0-5, The fraction 
of GR protons that enter into the CMZ, A, may be com¬ 
parable to /cmz- However, if the protons are not spheri¬ 
cally emitted and they are, for example, carried by strong 
outflows perpendicular to the disclike CMZ, the fraction 


may be much smaller. Moreover, the inner edge of the 
CMZ may not have contact with Sgr A*. Thus, we as¬ 
sume that A < 0.5, and A <C 0.5 is very likely. Observa¬ 
tions have shown that the mass of the SMBH is Mbh = 
4.3 X 10® Mq [^, and the current mass accretion rate is 
M = 4xlO“®MQyr“^ [^. Since the Eddington luminos¬ 
ity is given by TEdd = 1-26 x 10®® (MbhZ-Mq) ergs“^, the 
Eddington accretion rate is MEdd = 9.6 x 10“®M© yr“^. 
Thus, the normalized accretion rate is written as m = 
4.2 X 10“®. For these and the fiducial parameters, we 
obtain = 0.2 TeV from Eq. (P). We solve Eq. (P) 
from t = 0 to 10^ yr. There are no CRs at t = 0. The 
distribution of CRs has achieved a steady state at the 
end of the simulation because the diffusion time of CRs 
is only tdiff.c = 7icMz/(6'^) 1-6 x 10® yr at Ep ^ 1 TeV. 

The following results are those at to = 10^ yr. Since this 
time scale is much larger than tdiff.C) the energy spec¬ 
trum of the GR protons is almost uniform in the CMZ 
[i^ As long as the spectrum of the injected CRs 
does not vary significantly, our model does not expect 
substantial variation in the gamma-ray spectrum across 
the CMZ (see later), which is consistent with observa¬ 
tions. The distance to Sgr A* and the CMZ is assumed 
to be 8.5 kpc. 

The inverse of the cooling time of GR protons due to 
pp interactions is 


t 


-1 

pp 


'^CMTi^pp^^pp 5 


( 6 ) 


where ncMZ = PCMz/wp and Kpp[^ 0.5) is the proton 
inelasticity of the process. The total cross section of the 
process is given by 


CTpp = (34.3-f 1.88L-b0.25L2) 



(7) 


where Eth = 1.22 GeV is the threshold en ergy of produc¬ 
tion of TT® mesons and L = ln(Ep/l TeV) [^. Thus, the 
cooling time is tpp ^ 7 x 10® yr at Ep ^ 1 TeV, which 
is larger than tdiff,c- Since tpp > tdiff.c is satisfied in 
the energy of interest {Ep > 1 TeV), we do not need to 
include the cooling effect in Eq. ©. 

EigureP shows the gamma-ray flux and neutrino flux 
(per flavor) from the CMZ when rh = 4.2 x 10“® and 
A = 0.01 regardless of t. Other parameters are the fidu¬ 
cial ones. Eor comparison, we show the GeV and TeV 
gamma-ray fluxes obtained with Eermi and HESS obser¬ 


vations 


53l . l54j . The predicted gamma-ray flux is much 
smaller than the observations. The fraction of A ^ 0.5 
is needed in order that the flux is comparable to the ob¬ 
servations at E ^ 1 TeV. However, as is noted above, 
the value of A = 0.5 is probably too large for the actual 
CMZ. The gamma-ray and neutrino spectra are similar 
and they are not represented by a power-law, because 
they reflect the proton spectrum [Eq. ([2])]. 

Recent studies have indicated that the current activity 
of Sgr A* is exceptionally small, and that the average 
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FIG. 1. Predicted gamma-ray flux (dashed line) and neutrino 
flux (two-dot dashed line) from the CMZ when rh = 4.2 x 10“® 
and A = 0.01. Filled circles an d sq uares are the Fermi and 
HESS observations, respectively 1^. 



E (GeV) 

FIG. 2. Same as Fig.[T]but for rh = 0.001, and A = 5 x 10”^. 
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FIG. 3. Same as Fig. [T] but for rh = 0.001, A = 3 x 10 
C = 0.18, and rjcr — 6 x 10“®. 



E (GeV) 

FIG. 4. Same as Fig. \T\ but when rh is variable (see text). 
Thick dashed line is the total gamma-ray flux. Thin dotted 
line (out) is the gamma-ray flux from 0.9i?cMZ < r < Rcmz- 
Thin dot-dashed line (in) is that from 0 < r < 0.1 Rcmz- 


accretion rate more than 100 yrs ago might be much 
larger and it could be as much as 10^-10^ times the cur¬ 
rent one [H-ii. Thus, we calculate the gamma-ray and 


gardless of t. The drop of activity in the past ^100 yrs 
does not affect the results because the diffusion time of 
CRs is much larger than 100 yrs. Other parameters are 
the same as the fiducial ones. These give the typical en¬ 
ergy of ifp.eq = 3.4TeV from Eq. (P). Note that for given 
m and Mbh, any combinations of parameters (a, /3, 

and Rb.cc/Rs) that give the same Ep^eq give the same 
spectrum. Moreover, any combinations of A and ijcr that 
give the same Xrjcr give the same spectrum. Figure [2] 
shows the results of this model; the gamma-ray spectrum 
at E ^ 0.2-10 TeV well reproduces that obtained with 
HESS. If observations for E > 10 TeV become available 


in the future (e.g. Cherenkov Telescope Array; CTA ISSll . 
High Altitude Water Cherenkov detector; HAWC [56|). 
our model predicts a soft gamma-ray spectrum in that 
energy band. Since the designed sensitivities of CTA and 
HAWC at A ~ 10-100 TeV are better than that of HESS, 
the flux predicted in Fig. [2] could be easily detected. The 
gamma-ray image taken with HESS appears to coincide 
with the CMZ 5^, which supports this model. Since the 
apparent size of the CMZ is 3° x 0.5°, it can be well 
resolved by CTA with a resolution of ^ 1' [H^ . Detailed 
maps of gamma rays will reflect not only the distribution 
of molecular gas but that of CRs. The latter may reflect 
the history of Sgr A* activities, if the activities signifi¬ 
cantly change on the diffusion time scale of the CRs. 

From the current neutrino observations with IceCube, 
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the flux of ~ 3 X 10“® GeVcm”^ s“^ can be attributed to 
that from the Galactic center [l^. Since the predicted 
fluxes in Figs. [T] and [5] are smaller than that, they are 
consistent with the observations. However, as the statis¬ 
tics of neutrinos improve, we may detect an excess in 
the direction of the Galactic center in the future. In 
particular, observations with KMSNeT would be useful 
to detect the CMZ as a neutrino source if the flux is 
> 10“® GeV cm“^ s“^[53|- Our model predicts that the 
neutrino image should coincide with the gamma-ray im¬ 
age because both are the results of pp interactions. 

Since parameters for the RIAF have some uncertain¬ 
ties, we adopt another model that can reproduce the neu¬ 
trino flux at ~ 1 PeV obtained with IceGube [l^. We 
change the parameter for the turbulence to C = 0.18, 
and the acceleration efficiency of GRs in the RIAF to 
r^cr = 6 X 10“^. We take m = 0.001, which is the same 
as that in Fig. [21 and A = 3 x 10“^ so that the flux at 
E ^ 1 TeV is consistent with the HESS observations. 
Other parameters are the fiducial ones. For these pa¬ 
rameters, the typical energy is = 160 TeV (Eq. [T|), 
which is much larger than that of the model in Fig. [2] 
The results are shown in Fig. |3l If this is the case, a rel¬ 
atively hard gamma-ray spectrum would be observed by 
CTA at ^ 10-100 TeV, and could be discriminated from 
the spectrum in Fig. |21 The gamma rays at A < 1 TeV 
should have another origin. 

So far we have assumed that the accretion rate, m, 
is constant. Here, we discuss the effects of variable m. 
The x-ray light curve of Sgr A* in the past 500 yrs was 
derived in Ref. [1^. They showed that the x-ray lu¬ 
minosity is Lx ^ 10^® erg s“^ in the past 50-500 yrs, 
and then it dropped to the current value of Lx ^ 10^^- 
1035 gpg g-i^ qijjg x-ray luminosity before 500 yrs ago 
is less constrained. Upper limits are placed to down to 
about 8 X 10'*^°ergs“^ for several periods within the past 
4 X 10^ yrs 53. Before that, upper limits are ~ 10^^- 
10^2 grgs-i [Hall. The x-ray luminosity of a RIAF is 
proportional to ttt? 22|. Assuming that the x-ray lu¬ 
minosity follows the above observations and upper lim¬ 
its, and that to = 10^ yr is the current time, we set 
TO = 0.03 for 0 < t < to ~ 4 X 10"'^ yr, to = 0.01 for 
to — 4 X 10"* yr < t < to — 1 X 10^ yr, to = 0.001 for 
to — 1 X 10"^ yr < t < to — 50 yr, and to = 4.2 x 10“® for 
to — 50yr < t < to. Other parameters, including the fidu¬ 
cial values, are time-independent, except for A = 4 x 10“® 
and C = 0.025, which are chosen to be consistent with ob¬ 
servations at A ^ 0.2-10 TeV. 

Figure m shows the results. The gamma-ray flux from 
the outer region originates in CRs injected in earlier 
times. Since to decreases as time advances, the typical 
energy of the CRs, Ap^eq, in the CMZ should decrease 
from the outer region to the inner region [Eq. ([T|)]. We 
obtain Ep^eq = 2.3 TeV, when to = 0.03. However, the 
shape of gamma-ray spectrum from 0 < r < 0.1 i?cMZ is 
not much different than that from 0.9i?cMZ < r < i?cMZ- 


The peak gamma-ray energy of the former is only a fac¬ 
tor of 2 smaller than the latter. Most of the gamma-ray 
flux from the CMZ is associated with CRs injected when 
TO = 0.03, because they are injected during most of the 
past diffusion time (to—^diff.c <t< tg—4x10'^ yr), where 
tdiff.c = 4 ?cmz/( 5«) 1-6 X 10® yr (at Up — 1 TeV). 

Those CRs prevail in the CMZ even at present. CRs 
of U < 1 TeV injected when to < 0.001 are located 
at r ^ 0.25 i?cMZ- However, their short injection time 
scale compared with tdiff.c and the smaller injection rate 
(Tp,tot oc to) make their contribution to the gamma-ray 
flux smaller. In other words, the contribution of CRs is 
represented in the form of f XLp^totdt integrated for the 
past diffusion time. Moreover, since the higher-energy 
CRs (^ 1 TeV) have shorter diffusion times, they es¬ 
cape faster from the CMZ than lower-energy CRs. This 
also makes the gamma-ray spectrum from the older CRs 
softer. Of course, if the typical CR energy, Up^eq, sig¬ 
nificantly varies in the past, while to does not much 
decrease [see (Eq. [1])], the gamma-ray spectrum can 
change across the CMZ. Note that the gamma-ray lu¬ 
minosity of the RIAF is proportional to and it is 
^ 0.001 Lp,tot ~ lO^^erg s“^ when to = 0.001 [l^. Al¬ 
though this is larger than the current total gamma-ray 
luminosity of the CMZ at ^ 1 TeV (^ 4 x 10®^ergs“^), 
the gamma-rays from the RIAF cannot be observed at 
present. This is because the gamma-ray luminosity of 
the RIAF almost immediately changes with to, owing to 
the short (< 1 yr) diffusion or escape time of the CRs in 
the RIAF 0. 


V. GAMMA RAYS FROM THE CMZ AROUND 
CENTAURUS A 


Since some LLAGNs other than Sgr A* also have their 
own CMZs, gamma rays and neutrinos may be created 
there. As for neutrinos from their RIAFs, LLAGNs with 
TO ^ 0.01-0.1 most contribute to the neutrino flux on the 
Earth, because the neutrino production rate in a RIAF is 
sensitive to to and RIAFs are realized when to ^ 0.01-0.1 
0 . If a fraction A 10 ® of the CR protons acceler¬ 
ated in those RIAFs enter their CMZs as is the case of 
Sgr A*, the production rate of neutrinos in the CMZs is 
smaller than that in the RIAFs. Thus the contribution 
of the former to the overall neutrino flux on the Earth is 
expected to be smaller than the latter. However, if the 
RIAF in a nearby galaxy is well covered with massive 
molecular gas or a CMZ and to is relatively large, the 
gamma rays from the CMZ may still be detectable as an 
individual source. 

In Fig. |5l we show as an example the gamma-ray flux 
from Centaurus A (Cen A), which is a nearby radio 
galaxy and for which the origin of the gamma rays is 
under debate 62, The distance to Cen A is assumed 
to be 3.84 Mpc. We do not include the absorption of 
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FIG. 5. Predicated gamma-ray flux (dashed line) and neu¬ 
trino flux (two-dot dashed line) from Cen A. Parameters are 
shown in the text. Filled circles and squares are the Fermi 
and HESS observations, respectively 


gamma rays and neutrinos are created there through pp 
interactions. We showed that our model can explain the 
gamma rays observed by HESS at E ~ 0.2-10 TeV, if 
the accretion rate on Sgr A* was ~ 10^ times larger 
in the past than it is today as indicated by previous 
studies, and if the typical energy of the CR protons is 
~ TeV. In the near future, CTA could observe gamma 
rays at ~ 10-100 TeV, which could be used to estimate 
the typical energy of the CR protons more precisely. The 
gamma-ray emission from some nearby galaxies could be 
attributed to this mechanism if their LLAGNs are sur¬ 
rounded by molecular gas. Future comparison between 
observed gamma-ray and neutrino spectra and images 
would be useful to confirm this model. 

This work is partly supported by the challenge sup¬ 
port program of Osaka University and by JSPS KAK- 
ENHI Grant Number 15K05080 (Y.F). This work is also 
partly supported by Grant-in-Aid for JSPS Fellowships 
No. 251784 (S.S.K.). K.M. acknowledges Institute for 
Advanced Study for continuous support. 


the gamma rays. The radius, thickness, and mass of 
the GMZ are i?cMZ,obs =195 pc, TJcMZ.obs =195 pc, and 
Mcmz = 8.4 X 10^ Mq, respectively [H^. The mass of the 
SMBH is Mbh = 5 X 10^ M© [6i|. We choose to = 0.01, 
A = 0.02, and ( = 0.03 in order to reproduce the HESS 
results. Other parameters are the same as the fiducial 
ones. These give the typical energy of Ep^eq = 7.9 TeV 
from Eq. ©■ Cen A has a prominent cold gas disc [o^ 
and, thus, the effective covering factor A may be larger 
than that of Sgr A*. Figure [S] shows that our model 
can reproduce the HESS observations at A > 1 TeV, al¬ 
though another component is required at E 1 TeV. We 
note that the actual gamma-ray flux from the CMZ could 
be smaller, if that from the CRs accelerated by other 
mechanisms (e. g. acceleration in the electronic field at 
the base of jets |67|) cannot be ignored. In fact, jets have 
been observed in Cen A 681. 


VI. SUMMARY 

We have shown that TeV gamma rays from the 
Galactic center can be used to test a model in which 
low-luminosity active galactic nuclei (LLAGNs) are the 
source of neutrinos observed with IceGube. In this model, 
protons are accelerated in the radiatively inefficient ac¬ 
cretion flows (RIAFs) in the LLAGNs, and neutrinos are 
created through pp and pj interactions in the flows. Since 
Sgr A* at the center of the Galaxy is an LLAGN, we ex¬ 
pect that the protons are being accelerated in Sgr A* and 
injected into the interstellar space. 

In this study, we found that the central molecular zone 
(GMZ) surrounding Sgr A* works as an effective target 
of the high-energy protons escaped from the RIAF and 
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